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2 Abstract 

r-^ An exact description is provided of a flat fluid membrane locally deformed by external normal forces bringing 

^ two points on the surface together. The conformal symmetry of the two-dimensional bending energy is used 

c/) to identify the distribution of energy as well as the stress established in the membrane. Logarithmic curvature 



singularities in the surface geometry at the points of contact signal the presence of external forces. Whereas 
the total energy is independent of the distance S between the two points, as measured along the surface, the 
magnitude of the forces is shown to vary inversely with S. The geometry behaves near each of the singularities as 
a biharmonic monopole, in the region between them as a surface of constant mean curvature, and asymptotically 



I 



o 



as a biharmonic quadrupole. Radial tension is accompanied by lateral compression, both near the poles and 



O asymptotically, with a crossover from tension to compression displayed between the two singularities. 
> PACS: 87.10.Pq, 87.16. D-, 02.40. Hw 

cl Introduction 

ON 

^ Striking feature of cellular membranes in biology is that they behave as fluid membranes on mesoscopic length 
r^ales. Their mechanical properties are captured by the geometrical degrees of freedom of the membrane surface, 
^hich is described with extraordinary accuracy by the Canham-Helfrich bending energy, quadratic in the curvature 

labile it is simple to write down the energy, the non-linear field theory associated with it is not simple. Thus, unless 
the membrane lends itself to a description in terms of a height function which varies slowly above some simple 
geometry, it will tend to be necessary to examine its behavior numerically. 

On occasion, the symmetries of the theory allow one to circumvent this necessity. In this paper, a stratagem 
exploiting these symmetries will be used to provide an exact analytical description of a local deformation that is 
established in a flat membrane when two neighboring points are brought into contact. While this state describes 
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an idealization from a physical point of view, it is nonetheless a genuine non-perturbative feature of the mesoscopic 
description of fluid membranes, providing an unexpected access to the non-linear response of fluid membranes to 
external forces. 

Consider two points a fixed distance apart on an initially flat membrane. There is one obvious way to bring these 
two points together: apply tangential forces. In the geometric approximation there is no penalty associated with 
surface shear and thus the two points can be brought together this way without deforming the membrane. It is also 
possible, however, to apply a force normal to the surface in such a way that the points remain separated on the 
surface but are brought together in space. A cartoon illustrating this process is shown in Fig. [T] 




t t 



Figure 1: Sequence of normal deformations bringing two points on a flat membrane together. 
Two views of the end point of this process are illustrated in Figs. [2^ and |2]d. 

Even though the membrane undergoes large deformations in a neighborhood of the two points under the influence 
of the force, an exact description of this state exists, a happy accident of the conformal invariance of the two- 
dimensional bending energy ([21 El H])- In contrast, neither the initial state in the sequence-involving two external 
vertical forces-nor any of the intermediate state terminating in contact-with horizonal forces-lend themselves to an 
analogous treatment. This is because conformal invariance is broken in all but the final element in the sequence; 
for if the points are maintained a finite distance apart this distance introduces a scale that competes with the scale 
established by their separation on the surface. 

Axially symmetric membrane states with two poles brought together were constructed a few years ago [F, "G]. This 
was done by exploiting the fact that any state related by inversion in a sphere to an equilibrium state will also be an 
equilibrium state. In particular, states described by inversion of a catenoid-an axially symmetric infinite cylindrical 
surface growing exponentially at its two ends-are also equilibria. Under inversion, the topology of the catenoid 
will be changed as its two ends get compactified to a single point. The compact axially symmetric vesicle states 
described in [5l|6] result when the center of inversion lies within the catenoid. By positioning this point instead on 
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(a) top view: touching fingers (b) bottom view: membrane pinch 



Figure 2: Geometry of localized defect on sheet with two flat directions 

the surface of the catenoid its neighborhood get inflated into a large flat region. A very different beast, the localized 
dipolar defect, illustrated in Fig. [2] sits at the center of this flat geometry. 

From one point of view, the defect on the membrane resembles a flat sheet with two fingers of membrane-touching 
at a point- pushed up from below, a valley forming beneath them, as illustrated in Fig. [2^. A complementary 
description is also evident in terms of a flat sheet that has been pinched into self-contact; a ridge forms above the 
two points that have been brought into contact, illustrated in Fig. |2)d, not unlike the wrinkle that is formed in skin 
that has been pinched together between two fingers. One description is the other viewed from below. Note that 
a well defined tangent plane exists at the two points of contact. This is consistent with a logarithmic singularity 
in the curvature at these two points. These singularities originate in the compactification of the exponential ends 
IS Ej. The asymptotic behavior of the minimal surface completely determines the geometry at the center of the 
defect. The conformal symmetry of the bending energy thus provides a connection between the weak-field harmonic 
behavior of the catenoid at infinity, described accurately by perturbation theory, and the non-linear behavior of the 
defect at its center. This kind of correspondence is a general feature of field theories with an underlying conformal 
symmetry. 

The focus here, for simplicity, will be on symmetric states, arising when the point of inversion lies on the waist of 
the catenoid. Such states are characterized uniquely by the geodesic distance separating the two points along the 
surface. 

While the energy density of these states is localized around the defect center, the total energy has the unusual 
property that it is independent of the distance between the points. 

The normal force tying the points together provides a source of stress in the membranej^ If these forces were removed, 
the sheet would collapse into a featureless planar state, the only unconstrained equilibrium geometry with this 
topology. The stress finds itself concentrated in the neighborhood of the two points: it is also completely captured 
by the membrane geometry. The curvature singularities in the defect geometry are associated with these localized 
external forces, and manifest themselves as divergences in the stress. The strength of these forces is registered by 
an appropriate integral of the stress along a contour encircling one of the points. As will be demonstrated it is also 

^The catenoid generating this geometry is itself a stress-free equilibrium state of a fluid membrane, albeit one that is never realized 
physically. This is not a problem here for its only role will be to generate the state of physical interest. 
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possible to evaluate them exactly by deforming the contour to exploit the symmetry of the problem. They are found 
to vary inversely with the distance S between the two points measured along the surface. Conformal invariance 
also permits an exact analytic description of the inhomogeneous and anisotropic distribution of stress established 
in the membrane in the fully non-linear theory. This is the first instance-to our knowledge-of such a description in 
any non-trivial fluid membrane geometry. Three qualitatively different regimes are identified: the neighborhood of 
the singularities within the fingers, the valley that extends under them, and the asymptotic region. While a global 
description of the surface in terms of a height function does not exist, in each of these three regions it is possible to 
describe the surface in terms of its height above an appropriate plane, and to compare the exact stress tensor with 
its counterpart in the small gradient biharmonic approximation. In particular, it will be shown that, whereas near 
the singularities the geometry behaves as a biharmonic monopole, in the valley it behaves as a surface of constant 
mean curvature, and asymptotically it behaves as a biharmonic quadrupole. Characteristic distributions of stress are 
found to be associated with each region; both near the poles as well as asymptotically, radial tension is accompanied 
by an equal lateral compression. A crossovers from tension to compression is encountered in the region between the 
poles. 



2 Geometry of the Defect 

Consider a catenoid, with unit neck radius, aligned along the Z axisj^ Its polar radius p at a given value of Z is 
given by III [H] p{Z) = coshZ, so that the surface is described parametrically in the form: 

Xo{Z,ip) = {p{Z)cosip,p{Z)smip,Z) . (1) 



X 




(a) Catenoid (b) Defect 



Figure 3: (a) Surface curves on the catenoid obtained from its intersection with relevant planes: catenaries on the plane 
Y = with if = {Tq), vr (F^) are represented by red and green curves respectively, the circle Fn on the plane Z = in 
black, and curves on the plane X = 1 {jf) in blue and yellow. The center of inversion is indicated by the point i. (b) The 
image under inversion in a unit sphere centered at i. 



^Three-dimensional Euclidean space will be described by Cartesian coordinates X,Y,Z, with basis vectors i — (1,0,0), j — (0, 1,0) 
and k — (0, 0, 1). Polar coordinates p and cp are adapted to the X — Y plane. 
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The surface obtained by inversion of this catenoid in a sphere of radius R, centered at the point i on its neck, will 
be described by the new position vector X = i?^X/|Xp, where X = Xq — i, [7J, given explicitly byn 



~ R? 

X = ^^^-^^^^^^ (pcos(/? - l,psin(/?,Z) , P^{p,ip) = -2pcosip + l. (2) 
R^ has the dimensions of length. This is the defect geometry illustrated in Figs. [2]and|3| 

The two planar ends of the catenoid, represented by the neighborhood of the points at infinity Z ±00 get 
compactified to a single point, the origin of the Euclidean plane, O. 

This geometry possesses F = and Z = as mirror planes: the former is what remains of the axial symmetry; 
the latter is the unbroken mirror symmetry of the catenoid. These two symmetries are captured respectively by 

-H- —Lf and Z ^ —Z. 

To describe the defect geometry in greater detail, it is useful to trace the image under inversion of various surface 
curves on its two mirror planes: 

1. Consider the image of the two catenaries describing the intersection of the catenoid with the plane F = 
(given by y9 = 0, vr). 

The catenary with 93 = 0, Fq, passing through the center of inversion i, represented by the red curve in 
Fig|3|a), is mapped to two semi-infinite curves Fq and Fq (with Z < and Z > respectively) 

f^ = -—^-—Ap-ho,±\z\), (3) 

Z2 + (p - 1) 

both of which terminate at O and represented in red in FigjSfb). They both increase to an asymptotic value 

X = R^/2. 

The second catenary with = vr, F^r, maps to the finite curve F whose two ends also terminate at O. They 
are represented by the green curves in Fig|3][a) and (b) respectively. The point (—1, 0, 0) on the catenary maps 
to the point (—1/2,0,0) on this curve. Together F^ = Fg U F U Fq form a single curve with a continuous 
tangent vector, making contact with itself at O. 

2. The circular waist of the catenoid Fy resulting from its intersection with the plane Z = (illustrated with 
a black curve in Figj3][a)) passes through the center of inversion, and thus under inversion opens into the 
straight line f y (black line in Fig|3][b)) parametrized by 

f|| = ^(-l,cot(^/2),0), (4) 

which intersects F orthogonally. 

3. The curves 7=*" resulting from the intersection of the catenoid with the plane X = 1, given by taiup = ± sinh Z, 
get mapped to the curves 



^This surface is also parametrized by Z and (p. 
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which are the intersection of the defect surface with the plane X = 0. These two curves are represented 
in yellow and blue in FigsjSj^a) and (b). This plane represents the midplane of the asymptotic geometry, 
described in detail in section |2.2.3[ 



The curves t± and Fn are geodesies on the surface. 

The defect geometry is characterized by a single length scale, S, the length of the geodesic T connecting the two 
points in contact. This length is related to the radius of inversion R. To establish this relationship, note that the 
line element on the catenoid is given by 

ds^ = p\Z) [dZ^ + d^^) , (6) 

thus it is conformally flat, with conformal factor given by the polar radius p. Its counterpart on the surface, ds^ , is 
given by 

ds^ = ds^ : (7) 

thus the length of V is given by 

oo 

S = 2R-j z^^l^^y dZ ^ 1.206R^ . (8) 



A defect geometry with a given S is obtained from a unit catenoid by choosing R appropriately. For a given value 
of S there is a unique symmetric geometry. 



2.1 Curvatures 

The curvature of the defect will play a role not only in determining its energy but also in determining the stresses. 
Singularities in the curvature signal external forces. The curvature is characterized by two scalars, its extremal values, 
as well as the directions-mutually orthogonal-along which they occur. The intersections of the defect surface with 
the mirror planes are both directions of curvature as well as geodesic. 

On the catenoid, the maximum curvature is achieved along the parallel directions (circles of constant Z in Fig. 
|4][a), and the minimum curvature along meridians (catenaries in Fig. |4][a). They are given by = = —C±, 
confirming that the catenoid is a minimal surface, i.e. K = C\\ + C± = 0. Its Gaussian curvature K-g = C\\C± is 
given by 

)Ca = -l/p'. (9) 

The principal directions (along the parallel circles and along the meridians) are aligned with the adapted basis 
vectors = d^X. and = S^X, 

V|| = (— sinyj, cos V9, 0) ; V_l = 1/p (sinh Z cos v?, sinh Z sinyj, 1) . (lOs) 

In the appendix it is shown that the counterparts on the defect of these two vectors fields are given by an appropriate 
reflection: 

V, = RV,, 2=11, ± (11) 
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(a) Catenoid 



(b) Defect top view 



(c) Defect bottom view 



Figure 4: Principal curves and directions on the catenoid and on the defect, (a) On the catenoid, the meridians are catenaries 
(dashed lines) with tangent vector represented with blue arrows; the parallels (solid circles), with tangent vector Vy 
are represented by red arrows (they are oriented anticlockwise with respect to the symmetry axis). The normal vectors are 
represented by green arrows, (b) Top and (c) bottom views of principal curves on the defect. The catenaries get mapped 
to closed curves which make self-contacts at the poles, represented by dashed lines, whereas the parallel circles transform 
to rotated circles surrounding the poles, indicated by solid lines. Red, blue and green arrows represent the mapping of the 
principal directions and the normal vector onto the defect. 

where R = 1 — 2X (g) X, represents a reflection in the plane passing through the origin and orthogonal to X, and 
X denotes the corresponding unit vector. 

The integral curves of these two vector fields on the defect are represented in Figs. |4][b) and |4](c). In particular, 
observe that the parallel circles on the catenoid get mapped to two nested families of circles, each of which encloses 
one of the the poles; these circles degenerate into the straight line Tt, on the mirror plane, Z = 0. The meridians 
get mapped to a single family of nested closed curves, each of which passes through the point O where self contact 
is made. 

Three regions displaying features of particular interest can be identified on the defect; the neighborhood of the 
poles, the ridge/valley geometry connecting the poles, and the asymptotic region of the surface. Although the 
global defect geometry does not lend itself to a Monge parametrization involving a unique height function above a 
plane, it is possible to describe each of these regions separately in terms of a height function above an appropriate 
base plane. 

2.2 Height function representations 
2.2.1 Near region 

The region Z <^ 1 originates in the compactification of the exponential growing ends of the catenoid with Z ^ 1 
under inversion in a sphere centered on the point i. This region on the catenoid can be approximated by the height 
function hi = ln2p (p 3> 1) above the X — Y plane, {X = pcosLp and Y = psirup). Its counterpart on the defect 
under inversion in a sphere of radius R then has the height function represention 

hi = -R'p'\nU) , p=-, ^ = ^, (12) 
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above the plane Z = { X = i?^p cos and Y = R^psincj)). The two curvatures diverge logarithmically at the 
origin as p — )■ 

Cll,.^-^J^4+l±l) . (13) 



However their difference remains finite 



Cll-C^^-^^. (14) 



2.2.2 Ridge/Valley region 

The ridge/valley between the poles originates in the neighborhood of the point — i on the catenoid diametrically 
opposite the center of inversion. The latter is described by the saddle-shaped height function above the plane 
X = —1: h2 = — cos 20, where r = \/Y^~+~Z^, and (p = arctanF/Z. Its inversion in a sphere of radius R 
centered at i is a parabolic cylinder described by a height function above the plane X = —R^/2, {Y = i?^/4r sin0, 
Z = R^/A7cos4)) 

~ R^ ~ ~ 

/i2 = — cos^ , r = r, (f) = (f) , (15) 

or h2 = 4ZyR^ 

2.2.3 Asymptotic region 

The infinite region remote from the origin originates in the neighborhood of the center of inversion i on the catenoid. 
The latter is described by the height function = r^/2cos20. Inversion of this region in a sphere of radius R 
centered at i can also be described by the height function above the the asymptotic plane X = {Y = _R^rsin0 
and Z = i?^f COS0) 

~ R^ ~ 1 ~ 

/l3=— cos20, f=-, = 0, (16) 

2 r 

which is independent of the radial distance r on the base plane. Thus in the asymptotic region, with r 3> 1, the 
defect interpolates between the two orthogonal lines, X = R^/2, Y = and X = —R'^/2,Z = 0; its height set 
by the scale S (see Eq. (|8])). This geometry is illustrated in Fig. |5] Unlike the exact defect, the geometry is 
symmetrical (modulo a rotation by 7r/2) with respect to the plane X = 0. In the small gradient approximation, the 



Figure 5: The surface given by the biharmonic height function h = 1/2 cos 20 provides an accurate description of the defect 
in the asymptotic region. 
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bending energy is proportional to 

H^IJ d^{V'h)\ (17) 

where dx is the element of area on the plane and is the corresponding Laplacian. This is minimized when h is 
biharmonic on the plane, or (V^)^/i = 0. Note that all three height functions are biharmonic on their respective 
planes, consistent with the fact that each of the three surface geometries locally minimizes bending energy. The 
catenoid from which they stem is a minimal surface minimizing locally the surface areaj^ 

In cylindrical coordinates, h = cos20-describing the asymptotic region-is the only biharmonic function independent 
of r. It belongs, with A; = 2, to the family of biharmonic functions, h = cos kcp/r'^^'^ with a A;-fold dihedral symmetry, 
obtained by the inversion in a sphere centered at the origin of the harmonic height functions h^'^'> = r'^ coskcf), k >2, 
representing k-th order saddles (a monkey saddle is represented by A; = 3). 



3 Defect Energy 

In general, the bending energy of a surface is given by 

H[X] = t,H,[X\ + RHGBm, (19) 

where 

Hi[X.] = ^ JdAK\ HcBm = JdAKc, (20) 

in which the mean and Gaussian curvatures are given by the trace and determinant of the curvature tensor, K = 
Ci + C2 and Kg = C1C2. The microscopic parameters k and R are respectively the bending and saddle splay 
modulus, with units of energy. The catenoid is a minimal surface so that Hi = 0. Hgb is the Gauss-Bonnet 
topological invariant. 

While the bending energy is not itself invariant under inversion, it is modulo a contribution proportional to Hgb- 
This is because Hi can be cast as a linear combination of Hgb and the conformally invariant Willmore energy 
H\y = 1/2 j dA (C|| — C_l)^. As a consequence, equilibrium states are mapped to equilibrium states under inversion. 
The defect geometry constructed in section |2] is thus also an equilibrium state. 

It is simple to show that the bending energy of the inverted surface is given by (see also [5]) 

H = k{Hi + 2{Hgb - Hgb)) + kHgb. (21) 

Whereas Hi vanishes for a catenoid it does not on its inverted counterpart. Thus the total bending energy of a 
surface obtained by inversion in a sphere of a minimal surface is proportional to the difference of two topological 
contributions. 



In general the topology will change under inversion if the surface is infinite or the point of inversion lies on the 
surface. The catenoid has the topology of a punctured plane with Gauss-Bonnet energy Hgb = — 47r; the defect 

^The surface area excess is given in the small gradient approximation by 

A^^ldx(Wh)\ (18) 
so that minimal surfaces satisfy the Laplace equation locally. 
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topology is planar, with Hgb = (Appendix A). Thus the bending energy of the defect is given hy H = 8ttk, 
independent of R. 

In particular, H is independent of S. Although the total energy is constant, the distribution of energy in the 
membrane 



H 



-kK' + RK, 



G ■ 



(22) 



is concentrated in the vicinity of the defect center; it also depends on k. Note first that the mean curvature of the 



defect is proportional to the support function on the catenoid (see Appendix B) 

4 



K 



Xn 



Its (scaled) value, as well as K"^ /2, are illustrated in Figs. |6j(a) and [6 
cosh Z — Z tanh Z , indicated in black in FigjHJa) and[6jb). 



X ■ n = 1 — cos (fi sechZ — Z tanh Z . (23) 

It vanishes along the curves cosy? = 



Likewise, the Gaussian curvature on the defect is given by 

1 



G 



i?4 



4 (X ■ n)' - (Z^ + P(p, 0)2) ' sech^Z 



(24) 



Kg is positive in a neighborhood of the two poles that extends on its dorsal side to infinity, a region one would 
expect to be under tension; it vanishes where 4 cos (p = sechZ (Z^ + 1) + 3 cosh Z — 2Z sinh Z, and along the axis 



of symmetry Z = 0, indicated in black in Fig. J6|c). In section [2.2.1 it was seen that the two principal curvatures 
diverge logarithmically near the poles (see Eq.([13j), so that the bending energies do also: K'^,Kg ~ 41n2(p/2)/i?^. 
The corresponding densities, however, not only remain finite but vanish at the poles. Asymptotically, T-L decays as 



1/r , as is easily confirmed using Eq. (16) 




(a) tanh if 



(b) tanhl/2i4:' 



(c) tanh/Cc 



Figure 6: (a) K, (b) K'^/2, (c) Kq- These functions vanish along black lines. 
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4 Stresses and transmitted forces on the defect 



The stress tensor on a membrane, described by the energy H (19), is given by f" = + Z'' n (TUl E] (see also 
nil [121 H]), where 

P'' = K (^K"^ -^K g"^^ , r = - , (25) 

where e^, a = 1,2 are the two tangent vectors adapted to the parametrization, n is the normal vector, and gab and 
Kab are the the metric and extrinsic curvature tensors, defined in the appendix. Under a local deformation of the 
equilibrium surface, X — t- X + 5X, its energy changes by an amount 

5H = - j t^AP-Va^X. (26) 

This stress depends only on the local geometry. There is no local stress associated with the Gauss-Bonnet bending 
energy. The tangential stress is quadratic in Kab and thus assumes its maximum and minimum values along the 
directions of curvature. In a source free region, P is conserved in equilibrium so that [lOj 

V,r = 0, (27) 

where Va is the covariant derivative compatible with the metric tensor. 

Consider a curve 7 dividing the surface in two regions, such as one of those illustrated in Fig. [7] The force per 
unit-length on one side of the surface (darker region in Fig. [T]), exerted by the other side (lighter region in Fig. [T]) 
across a line element of 7, is given in terms of the stress tensor by = — /af", where /„ = gabl^ are the covariant 
components of the vector 1, the outward normal to 7 tangent to the surface, that is, directed towards the region 
exerting the forcej^ If 7 is closed, the total force is given by the line integral 



dsi±, (28) 



where s is arc-length along 7. If 7 is contractible to a point or if it does not encloses a region with a singularity 
related to a source of stress, this integral vanishes. 



4.1 Stress distribution in the defect: quadratic approximation 

Before attempting to describe the fully non-linear distribution of stress in the defect, it is useful to first provide 
an approximate description of this distribution in each of the three regions described in section |2.2| in terms of a 
height function. It is simple to check that gradients remain small in the regions of interest, i.e. |V/i| ^ 1 where 
V represents the gradient on the base plane, so that a quadratic approximation to the stress tensor in terms of 
gradients and Hessians of the height function is legitimate. 

In the quadratic approximation, the tangent and normal vectors to the surface can be decomposed with respect to 
an orthonormal basis adopted to a Cartesian description of the appropriate base plane {E^,Ej^, k}: = Ej + Sj/ik 

^If t is the unit tangent vector and n the surface normal along the loop, then 1 = t x n and {t, n, 1} form a right-handed trihedron 
(the Darboux frame). Since t and 1 are tangent vectors to E they can be expanded with respect to the adapted basis vectors, {ea} as 
t = Tea and 1 ^ l^e^. 

^To unburden the notation, set i? = 1 and rotate the geometry so that the base plane coincides with the plane Z = Q parametrized 
by the radial and azimuthal coordinates r = y/X'^ + Y"^ and (j) — arctan(y/X). 
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and n = —dihEi + The metric tensor is given by gij = 5ij + VihVjh. In this approximation, the extrinsic 
curvature tensors is given by minus the Hessian of the height function Kij = —ViVjh. 



The quadratic approximation to the stress tensor (25) is now given by P = T^-'Ej + A^'k (compare with ^5J), where 



= v'/i {v'V'h - ]^{V^h)5''^ - V\V^h) V'h , N' = V'V^h . (29) 
Both T*-' and A^* are conserved, or VjT*-^ = and VjA^* = 0, when the height function is biharmonic. These 



equations are equivalent to the conservation law (27) in the quadratic approximation 



Let T and L represent the unit vectors along the projections of the vectors t and 1 (introduced below Eq.(28)) onto 



the base plane, so that they form an orthonormal basis adapted to the curve F obtained by projection of the surface 



curve 7 onto this plane. According to Eq. (28), the force per unit length transmitted across the curve 7 on the 
surface is given, in this approximation, by = —LiP = —LiT^^Ej — LiN'^k. 

4.1.1 Poles: /ii = -r2ln(r/2), r < 1 

Note that V^/ii = — 4(ln(r/2) + 1), and (V^)^/ii = — 87r5(r). Thus, near the poles, the stress tensor on the the 
base plane has components 

T''" = -41n^, T'^^ = T^'^ = 0, T'^^ = -lln^; (30) 
with counterparts normal to the plane 

N'' = --, N'^ = 0. (31) 

r 

Construct a circle of constant Vc on the base plane, oriented such that T = E,^ = (— sin0, cos 0, 0) and L = 
— Ej. = (— COS0, — sin0, 0)|^ By projecting the stress tensor onto L, one identifies the force transmitted per unit 
length across this circle to be given by 

fx. = -4(^ln|E, + ik^ . (32) 

Thus, the external radial force, viewed from the base plane, can be decomposed into a sum of radial compression 
and a vertical force associated with a source pulling (down) on the pole at its center. This is a striking illustration 
of the shortcomings of the Monge representation: one would expect radial tension rather than compression. It 
should be remembered, however, that E^ is not the tangent to the surface, nor is k its normal. To resolve this 
discrepancy, it is necessary to expand the radial force with respect to a basis of tangent vectors adapted to the 



surface: = E^ — 2r ln(r/2)k and n = 2r ln(r/2)Er + k. One then recasts Eq.(32) as 



l_Lr = 4 ( In - ^ n ) , (33) 



^The normal points toward positive Z 



®Note that T^, unlike the genuine tangential stress f"'', is not symmetric. Consequently, V iTji ^ 0. In addition, its trace T\ is 
different from zero, whereas /"^ vanishes. 

^This gives the force exerted by the interior region with r < Tc, containing the pole, on the exterior region with r > Vc 
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which indicates that the outer region is subjected to radial tension plus a normal force pulling on the pole. This 
coincides with the result obtained below in the non-linear analysis, Sect. 4.2 



If ij_r is integrated along the circle one identifies the total external force acting on the pole to be — Syrk. In fact, 



this quadratic analysis provides the correct exact non-linear result (42) derived below. This is because the quadratic 
approximation becomes exact at the pole. 

Now consider a radial line on the base plane along direction E,. and with outward normal L = E,J^ The force per 
unit length transmitted across this line is given by 

f|l^ = 41n^E^, (34) 
which implies compression along the azimuthal direction. 

Near r = 0, the azimuthal compression is equal in magnitude to the radial tension, as measured along the surface. 
This result hints at a more general result: tension in one direction is always accompanied by compression in an 



orthogonal direction. As discussed in greater detail in section 4/2 this is a direct consequence of the vanishing 
trace of the tangential stress, i.e. f\ = 0. Surprisingly, this is not evident if tangents are replaced by their 
projections on the plane, where one registers compression in both directions, radial and azimuthal. As was seen, 
this apparent discrepancy is an artifact of the Monge parametrization; specifically the third term in the definition 



of T*^, Eq. (29), originating in the projection of the normal stress onto the base plane, contributes a trace to T'-': 
Tii = — Vj(V^) Vj/i. Furthermore Vi{'V'^h)'Vjh diverges near the pole as 1/r. Thus, whereas the distinction 
between n and k disappears near the pole, where the gradients of h are small (n — k ^ 2r lnr/2 e^)), T*-'(n — k) 
diverges logarithmically. This example illustrates how the Monge representation can yield counterintuitive results 
whenever curvature singularites are involved. One needs to be careful to distinguish between tangent vectors and 
their projections onto the plane. This kind of gauge artifact is sidestepped in the covariant non-linear analysis 



presented in section 4.2 



4.1.2 Valley: /i2 = 1/4X2 

The tangential stress is diagonal with constant components 

Txx = l/8, Txy = Tyx = 0, Tyy = -1/8. (35) 

Furthermore, because the Laplacian is constant, V'^h2 = 1/2, both normal components vanish, Nx = Ny = 0. 
The force transmitted across the valley (constant X, along direction E^ with L = E^) is given by i±x = — 1/8E^, 
so that it is under compression. Likewise the force along the valley (constant Y along direction E^,, with L = E^) 
is f||y = l/SEy so that it is under tension along its length. 

4.1.3 Asymptotic region: /is = 1/2 cos 20 

In the asymptotic region, the components of the stress tensor on the plane are {V^h^ = — 2/r2cos20) 

T'"- = -^cos^2^, T"<^ = -T'^'- = ^sin40, T"^* = 4(1 + sin^ 20) . (36) 
^°The force is that exerted from a region with a value of </> below a given value on that with a higher value. 
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T*^ is not diagonal. This reflects the fact that the principal directions do not coincide asymptotically with the radial 
and azimuthal directions. The corresponding components normal to the plane are 

Ar^ = 4cos20, A^'^=4sin20. (37) 

Projecting onto the azimuthal and radial directions, one finds the forces per unit length transmitted across curves 
with constant r (L = E^,) and (L = E^) are given respectively by 

1 - 2 - A 

fj_r = — T sin 4<^ E^ + — cos^ 20 Er r cos 20k, (38a) 

^4 ^ 1^4 lyO 

2 -1-4 

f|l^ = -— (1 + sin^ 20)E^+— sin40E^- — sin20k. (38b) 

The tangential forces decay as r"'', their normal counterparts as r^^ , confirming the localization of the stress in the 
neighborhood of the poles; even though the geometry itself is not asymptotically flat, the stresses that unpinning it 
vanish asymptotically. Notice that, whereas the membrane is everywhere under radial tension far from the defect, 
its magnitude |fxr| varies from a maximum along the directions of the external forces (0 = 0,7r) to a minimum in 
the orthogonal directions (0 = 7r/2,37r/2). 

4.2 Distribution of stress: exact results 

The distribution of stress in the defect is completely encoded in the catenoid geometry. It is straightforward to 



show, using results presented in Appendix B that (see also [9J) 



f " = (^) ' (r f + ^ [k^' - ^^^?"') w,) , (39) 

where [H [T7] 

w„ = X X (n X e«) . (40) 

Because is orthogonal to X, Rw^ = w^. Whereas the stress tensor P vanishes in a minimal surface with 
= 0, it will not vanish in its inverted counterpart: 

f'^ = 4^i^"'w,. (41) 

It can also be shown (see [9]) that if is conserved then f" is also, except at the points in contact at X = 0, 
represented by Z — )• ±oo. It is evident that the stress is concentrated in the neighborhood of these points. Its 



divergence indicates the presence of sources on the right hand side of Eq. (27). These sources are identified with 
the external forces enforcing the contact constraint. In the absence of these forces the equilibrium would collapse 
to a plane with a vanishing stress. 
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The presence of a component along in fj^,., as well as component along E,. in fp^ reflects the fact that the vectors E^jE^ and 



k do not coincide asymptotically with the physically significant vectors, Vj^, Vii, and n, described in section 2.1 
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4.2.1 External forces on poles 



The external force pulling the poles together is determined by evaluating the integral of the stress Eq. (28) along a 
loop enclosing each point [18JP An appropriate contour is one of the closed circular integral curves of V\\ (indicated 
using red arrows in FigQb) illustrated in fig|7][a)j^ The mirror symmetry in the plane Z = can now be exploited 





(b) 



Figure 7: Closed integral curves of V|| and V± 



to deform this circular contour into the straight line Fji running along Z = illustrated on Figs. 3](b) and 4][b). 
The contribution to the force from the contour at infinity vanishes due to the fast asymptotic decay of the stress 
tensor (see Eq. (38a)). The external force on the right hand pole {Z > 0) is now given 



2n 



dslz 



i?2 



d^p (1 - COS(y9) = k. 



(42) 



The detailed distribution of stress need not be known to determine this force. This is just as well, for the detailed 
distribution of stress is not simple. This result also validates the analysis of the linearized theory in section 4.1.1 
where a circular contour, collapsing to the point of application of the external force, was used to determine this 
force. 



Note that there is no net vertical displacement of the membrane. This is consistent with the fact that the forces 
acting on the membrane are horizontal. Raising the membrane with respect to the asymptotic plane X = would 
require the introduction of vertical forces which would then need to be balanced on the boundary of the sheet. 

■"^^The forces supporting the defect cannot be resolved by examining the asymptotic physics. A contour completely encircling the 
defect will register a vanishing total force. 

■"^^The appropriately oriented boundary has the Darboux basis t = signZ V|| and 1 = signZ Vj^. 

■"^^The outward normal to Fn traversed from oo (tp = 0) to — oo ((^ = 27r) is 1 = k. In particular the components of the conormal 
vector expressed in the tangent basis are = |Xp/i?^, |Xp = 2(1 — cosip) and l'^ = 0, so that the only component of its associated 
covector 1 is = l/l^ . Also the line-element along Fn is ds = i?^/|Xp(i(/3 and the only relevant component of the stress tensor is 
given by = 2\X\^/R^k. 
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Using Eq. (8j) which relates the inversion radius R to the length S between the poles, one obtains FS = 9.65397r. 
This force diverges as the two points merge along the surface. 



It was seen in section [3] that the defect energy is independent of the distance S. If the constraint on S is relaxed, 
the two points become free to migrate across the surface by the process of conformal diffusion [4]. In due course, 
they will-on occasion-approach each other along the surface with a corresponding increase in the force. Below 
some critical separation, the maximum normal tension supported by the membrane or by the agent tethering the 
two points will be exceeded and the defect will disintegrate, decaying into a flat membrane as its bending energy is 
dissipated. Consistent with this picture, one would also expect a defect of a given size to possess a finite lifetime. 
This lifetime will, however, depend on physical properties of fluid membranes that fall outside the scope of the 
geometric description provided in this paper. 



4.2.2 Stress distribution revisited 



Using Eq.(39), one finds for the force per unit length transmitted across a circular loop with tangent Vn is given in 



terms of quantities measured along the parent catenoid by 



4|X 



2 



f± = -Va±r = ^C'^X X V|| . (43) 

A similar identity is obtained for fy with ± replaced by ||. The stress is clearly localized near the poles, represented 

by |X| — )• oo. 

Although the exact expression for the stress tensor is simple to write down, its physical interpretation is not transpar- 
ent. This is facilitated by examining its projections onto the orthonormal basis adapted to its curvature, {V^, Vy , n} 
discussed in section [2711 One has 

i±:=h±y± + f±n, (44) 
where f_L± = -Vj.aVj.bf'' and f± = -Vj_a/" are given by 

f±± = ;^C'x|X|2X-n, (45a) 

f± = -|jC^|XpX-V^. (45b) 

The counterpart of f^^, f|||| satisfies fjm = — f_L_L, a direct consequence of the traceless character of the tangential 
stress tensor, which is itself a manifestation of the scale invariance of the (Willmore) energy. The off-diagonal 
projection fy^ = V||aVj_bf"'' vanishes in this frame^^ If the membrane is under tension in the direction V^, it 
will be under an equal and opposite compression in the orthogonal direction V||, and vice-versa. The tension and 
compression will also assume their maximum and minimum values along these directions, f^^ is represented in Fig. 
8][a) for the principal curves Vy represented in Fig. 7^3). It is proportional to the mean curvature of the defect 

geometry: f^^ = ~(|X|/p)^-ft'; thus tangential stresses vanish wherever K = (cf. Fig. [6]^a)). Tension and 
compression anticorrelate with the sign of the mean curvature. 

Near the poles the defect is under radial tension in response to the external force. Moving away from the poles, 
the dorsal side of the pole remains under radial tension; however, a transition from radial tension to compression 



■"^^There is, of course, no shear stress supported by a fluid membrane. A non-diagonal tensor will occur when the tensor is decomposed 
with respect to an inappropriate frame. 
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(f±± < 0) occurs upon descent into the valley that stretches between the poles|^ The radial tension is accompanied 
by lateral compression, and the lateral compression along the valley is accompanied by tension along its lengthj^ 
The membrane is under rapidly decaying radial tension everywhere far from the defect center. These results are 



consistent with those obtained in the quadratic approximation in section 4.1.3 



The projection of the normal stress f± is represented in Fig. 8|b). It is negative everywhere, except along Fy where 
it vanishes, so that the normal external force transmitted across these closed curves points in a direction opposite 
that of the surface normal (see Fig. |4](c)). 



18 



As illustrated 



For completeness, the counterpart /y normal to the principal curves is represented in Fig. 8|c) 
in Fig. yj^c)), /|| is concentrated within the valley region and positive everywhere (directed along the normal) to the 
surface except along T±, where it vanishes. 




Figure 8: tanh of stress distributions in defect. Black lines represent points where they vanish. 



5 Discussion 

In general, the non-linearity of the underlying geometrical theory makes it very difFicult to obtain analytical insight 
into the physical behavior of fluid membranes. The conformal invariance of the bending energy, however, provides 
a window into their behavior-albeit idealized-under the influence of external forces, a regime one would not have 
expected to admit an analytical description. 

This symmetry has been used to construct a model of a fluid membrane with a local defect at its center held together 
by external forces acting normally to the membrane. This model provides a template for exploring the behavior of 

■"^^This asymmetry in the distribution of stress under external forces has its counterpart in a familiar daily routine: pulling a T-shirt 
by the back of the collar over one's head; whereas tension is established in the shirt along one's back, it will be compressed along one's 
chest. 

■"^^The distribution of stress along the valley is also consistent with one's expectations for a surface of approximately constant positive 
mean curvature. 

■"^^The appropriate Darboux basis, indicated in Fig. 7[b), is given by t = — signF Vj^ and 1 = signy V|[. 
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fluid membranes when localized external forces are applied. The description is non-perturbative. The deformations 
of the planar geometry are large at the center of the defect and do not lend themselves to a global perturbative 
description. 

More abstractly, the conformal symmetry implies a duality between the weak far field description in one geometry- 
the catenoid-and the strong near field behavior in another-the defect. Whereas the catenoid is a stress-free state, 
its dual exhibits a localized distribution of stress, supported by external sources, with curvature singularities at its 
center. At a technical level, this duality establishes a link between harmonic functions which describe the local 
behavior of minimal surfaces and biharmonic counterparts which describe an unconstrained fluid membrane. 
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Appendix A Gauss-Bonnet Energy 

Here the transformation of the Gauss-Bonnet energy under inversion is examined. Specifically it will be shown that 
the change under inversion is a sum of boundary terms. To reduce the burden of notation, consider the inversion in 
the unit sphere. One has 

dA}CG= I dAlCG-2 I dA (1^) K + A I dA (|^) • (A.l) 



However, 



X 2 / /V X 



|2 



X-n\^ 1 



Thus 

' dAlCG= [ dAJCc- [ dAVHn \X\\ (A.3) 



In any axially symmetric geometry, the Gauss-Bonnet energy itself can be cast explicitly as a boundary term. To see 
this note that Ci = 0, and C2 = sinQ/R where R is the polar radius, is the turning angle, and dot represents 
differentiation with respect to £, the arc-length along the meridian [7]. Replacing £ by as variable of integration, 
one obtains Hgb as a difference of cosines 

Hgb = 271 diRe{i) = 27r (cos 02 - cos 0i) . (A.4) 

Jo ^ 

The Gaussian energy of a surface of revolution depends only on the boundary values 0i and 02 of 0. For a 
catenoid, Hgb = — 47r; whereas Hgb = 47r for its axially symmetric inverted counterpart of spherical topology with 
a well-defined tangent plane at its poles 
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■"^^The sphere of inversion is not centered on the surface, otherwise the resulting surface has a planar topology and in consequence 

Hgb - 0. 
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This discrepancy is consistent with Eq.(A.3). To see this note that boundary term on the right in this equation is 
a sum of two terms of the form _ _ 

2 I dJ^-^^ = -An , (A.5) 

one for each pole. 

Appendix B Transformation of the geometry under inversion 

Let be the two tangent vectors to the surface adapted to the parametrization by Z and Lf, The corresponding 
tangent vectors on the defect geometry are given by 

ea= Re,, (B.l) 

where the linear mapping R = 1 — 2X ® X, represents a reflection of surface points in the plane passing through 
the origin and orthogonal to X. X denotes the corresponding unit vector. 

The metric tensor induced on the surface are given by gat = ■ e^. The metric under inversion is related to gab by 

gab gab = gab ■ (B.2) 

Let n be the unit normal to the catenoid, and Kab = ■ d^n the extrinsic curvature tensor on this surface. The 
counterparts on the defect are given by fi = -RnQand 

Kab = (Kab - ■ "^dab^ ' (^-3) 

It involves only the two tensors, gab and Kab- A derivation is provided in reference |9]. Let C\\ and C± denote 
the principal curvatures and Vy = Vj"ea, = V^e, the corresponding directions on the catenoid, so that 

K%Vl' = QV^, i =11, ±. Then K%Vf' = CiV^, where 

The corresponding principal directions are the normalized images of Vy and on the transformed surface: 

V, = i/"S„ = RV,. (B.5) 

These directions are thus appropriate rotations of their counterparts on the original surface and their components 
transform as V," = \'X\^ / R'^VI' . 



As a consequence of the relation (B.4), the difference of the two principal curvatures under inversion is preserved. 



modulo a scaling by the distance function and a change of sign. 



IX 



2 



qi-c^± = -^(qi-C7j. (B.6) 

In particular, on an inverted minimal surface (with a vanishing mean curvature, K = 0), one has K = 4/i?^X-n. 



minus sign preserves the orientation of the adapted basis. 
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